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Abstract 


In this study, analysis of the mitochondrial and Y-chromosomal haplogroups 
was used to provide information about the phylogenetic groups of Tutankhamun's 
family members, and their presence among the reported contemporary Egyptian 
population data. The mitochondrial and Y-chromosomal DNA analysis of 
Tutankhamun's family confirms our previous data ofthe royal family pedigree, 
with multiple controls authenticating all results. The proposed sibling 
relationship between Tutankhamun's parents, KV55 (Akhenaten) and KV35YL, 
is further supported. The royal lineage is composed of the Y-chromosome 
haplogroup R1b and the mitochondrial haplogroup K. Population genetics 
point to a common origin at ca. 14. 000—28. 000 years before present locating 
to the Near East. 


Introduction 


The analysis of mitochondrial DNA (mtDNA) has been an approach of interest to study the history 
of populations' migration and phylogeny (Macaulay et al. 1999, 232—249; Cavalli-Sforza and Feldman 
2003, 266-275; González et al. 2006, 124; Torroni et al. 2006, 339—345; Baca and Molak 2008, 
39-61). Compared to nuclear DNA, the hypervariable regions of mtDNA accumulate mutations 
rapidly. Being non-recombining and maternally inherited, mtDNA has become the molecular marker 
of choice to investigate genetic flow, to study human evolution, and to construct distinct phylogenetic 
groups (Giles et al. 1980, 6715—6719; Richards et al. 1996, 185—203; Macaulay et al. 1999, 232-249; 
Cavalli-Sforza and Feldman 2003, 266—275; Pakendorf and Stoneking 2005, 165-183; González et 
al. 2006, 124; Torroni et al. 2006, 339—345). 
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However, mtDNA identifies only the female perspective of a population and is liable to genetic 
selection (Pakendorf and Stoneking 2005, 165—183). Therefore, it has been integrated and supported 
with the study of paternally inherited Y-chromosomal DNA. This latter has been used to establish a 
global Y-chromosome phylogenetic tree analogous to that of the mitochondrial one (Jobling and 
Tyler-Smith 2003, 598—612; Baca and Molak 2008, 39-61). 


One of the most renowned royal families in ancient Egypt was that of the late 18^ Dynasty (ca. 
1550-1295 B.C.); including the pharaoh Akhenaten and his successor, Tutankhamun. The identity 
ofthe mummies, the relationship between them, and the description of inherited genetic diseases 
and pathological hallmarks in this dynasty have not been fully elucidated. Also, it is still not 
known to what extent endogamy has contributed to the early demise of the young king Tutankhamun 
and the consequent ending of the 18" Dynasty. Genetic fingerprinting together with statistical 
analysis was used in our previous study, to reveal the paternal and maternal autosomal lineages 
of Tutankhamun's family (Hawass et al. 2010, 638—647). The inherited autosomal half-allele, 
the partial Y-chromosomal data, and the possible trio (mother, father and child) figuring of 
Tutankhamun's family warranted the assembly of the most possible pedigree of the king's direct 
lineages. Yuya and Thuya were confirmed as the great-grandparents of Tutankhamun (generation 
I). Amenhotep III and KV 35EL were identified as his grandparents (generation II), and the KV 
55 male and the younger lady in KV 35 (KV35YL) were found to be his sibling parents (generation 
HD. Consolidating these findings with the historical records, KV 55 could be identified as 
Akhenaten and KV 35EL as Queen Tiye (Hawass et al. 2010, 638—647; Habicht et al. 2016, 
8216-8231). 


To expand the knowledge of the royal pedigree, another study was conducted to investigate the family 
maternal and paternal lineages. The current study employed analysis of the mitochondrial and 
Y-chromosomal haplogroups to test our previous results, (Hawass et al. 2010, 638—647) and to learn 
more about the genealogy of Tutankhamun's family. 


Materials and Methods 
DNA extraction, authentication and contamination monitoring: 


Strict adherence to requirements of ancient DNA authentication, including laboratory design and 
practice standards, was adopted for proof of reliability and quality of the results (Richards et al. 1995, 
291—299; Roberts and Ingham 2008, 600—613; Hawass et al. 2010, 638-647; Hawass et al. 2012, 
e8268; Keller et al. 2012, 698; Khairat et al. 2013, 309—325). DNA sampling and extraction were 
done following the published protocols for ancient DNA (Scholz and Pusch 1997, 61—64; Hawass et 
al. 2010, 638-647; Hawass et al. 2012, e8268). Sampling, DNA extraction, amplification, and 
sequencing were carried out in distinct working areas. Surfaces and equipment were continuously 
decontaminated with DNase and bleach and exposed to UV irradiation. Procedural components were 
strictly controlled and tested for the presence of contaminants. Negative controls were processed 
alongside each procedure to detect possible contamination with exogenous DNA sources (Tully et 
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al. 2001, 83-91; Pusch and Bachmann 2004, 957—964; Gilbert et al. 2005, 541—544; Lamers et al. 
2009, 40—55; Haak et al. 2010, e1000536; Hawass et al. 2010, 638—647; Fulton 2012, 1—11; Khairat 
et al. 2013, 309—325). DNA fingerprinting for the lab members was carried out to exclude the presence 
of staff DNA contamination. Reproducible results were verified through bacterial cloning of amplicons 
(Pääbo et al. 2004, 645—679; Fernández et al. 2014, e1004401) and through the analyses of different 
biopsies per mummy, where 3—6 bone samples were taken from various body areas of each studied 
mummy (Hawass et al. 2010, 638—647). Recruitment of a dedicated laboratory separated from the 
original one, run by an independent research team, was used for replication purposes (Gilbert et al. 
2005, 541—544; Lamers et al. 2009, 40—55; Hawass et al. 2010, 638—647; Fulton 2012, 1—11; Hawass 
et al. 2012, e8268). 


Mummies and subjects: 


Ten mummies were investigated in this study (Table 1), including Tutankhamun, his identified parents, 
grandparents, and great grandparents. In addition to the analyzed 10 mummies, the published data 
of a father-son pair [Ramesses III and Unknown man E/Pentaware] (Hawass et al. 2012, e8268) were 
included for haplotype comparison. 


Y-STR profile analysis: 


The AmpFISTR Yfiler PCR amplification Kit (Applied Biosystems, Foster City, California) was used 
to analyze the Y-chromosome specific haplotype profiles. The kit included sixteen human Y-chromosome 
DNA markers with sizes ranging between 99 bp and 324 bp, and a mean amplification size of 200 
bp. All markers were tetranucleotide repetitive microsatellites with the exception of DYS392, DYS438 
and DYS448, which displayed, respectively, trinucleotide, pentanucleotide and hexanucleotide repeat 
motifs. DNA templates were amplified according to the manufacturer's protocol with some modifications 
that involved decreasing the concentration of primers by 40% and increasing the AmpliTaq DNA 
polymerase concentration by 50%. Detection of the amplified PCR products was carried out using 
the ABI Genetic Analyzer 3130, Data collection Software v3.0 and GeneMapper ID v3.2 (Applied 
Biosystems). 


The majority rule was employed for STR allele designation, where a predominant peak was generated 
at least 10 times [i.e. each exceeding the threshold of 50 relative fluorescence units] per mummy, 
using various bone samples from different body areas and monitored for slippage (Hawass et al. 2010, 
638—647). A final consensus profile of each investigated mummy was constructed using the reproducible 
data of the generated partial profiles (Cowen et al. 2011, 400—406). 


Amplification of the mitochondrial hypervariable region I, cloning and Sanger sequencing: 
Specific primers were designed to amplify part of the hypervariable region I (HVRI) encompassing 


positions 16056 to 16409 (Table 2). These mitochondrial segments were amplified in a 50 uL PCR 
mix containing 1 X AmpliTaq Gold® PCR Master Mix (Applied Biosystems), 20 pmol of each primer, 
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Name|Primer Pair Sequence Annealing Product 
(5'-3') Temperature | Size(bp) 

m10 | L-15931 ACCAGTCTTGTAAACCGGAG 59°C 250 
H-16141 ATGTACTACAGGTGGTCAAG 

m8 | L-16056 GAAGCAGATTTGGGTACCAC 59°C 220 
H-16232 (GCTTTGGAGTTGCAGTTGATGTGT 

m24 | L- 16059 AAGCAGATTTGGGTACCACCCA TDP*62-57-52 97 
H-16111 [GGTACCGTACAATATTCATGGTG (5,5,30) 

m3 | L-16118 [TACATTACTGCCAGCCACCAT 62°C 159 
H-16232 [GCTTTGGAGTTGCAGTTGATGTGT 

m11 | L-16118 TACATTACTGCCAGCCACCAT 55°C 225 
H-16302 |TGGCTTTATGTACTATGTAC 

m4 | L-16118 [TACATTACTGCCAGCCACCAT 59°C 140 
H-16217a [TGTGTGATAGTTGAGGGTTG 

m18 | L-16126 CTGCCAGCCACCATGAATATTG 60°C 185 
H-16266 |TAGGTTTGTTGGTATCCTAGTGG 

m23 | L-16129 (GCCAGCCACCATGAATATTGTAC TDP62-57-52 134 
H- 16217b TGATGTGTGATAGTTGAGGGTTG (5,5,30) 

m26 | L-16129 GCCAGCCACCATGAATATTGTAC 64°C 150 
H-16232 [(GCTTTGGAGTTGCAGTTGATGTGT 

m12 | L-16171 |CCACCTGTAGTACATAAAAACCCA 63°C 109 
H-16232 |GCTTTGGAGTTGCAGTTGATGTGT 

m6 | L-16210 CCCCATGCTTACAAGCAAGT 529C 132 
H-16302 [TGGCTTTATGTACTATGTAC 

m5 | L-16210 CCCCATGCTTACAAGCAAGT 62°C 178 
H-16347 |ATGGGGACGAGAAGGGATTTG 

m17 | L-16211 |CTCCCCATGCTTACAAGCAGTA 609C 100 
H-16266 JAGGTTTGTTGGTATCCTAGTGG 

m22 | L- 16226 GCAAGTACAGCAATCAACCCTC TDP62-57-52 109 
H- 16290 GTACTATGTACTGTTAAGGGTGG (5,5,25) 

m7 | L-16288 CACTAGGATACCAACAAACC 62°C 161 
H-16409a GCGGGATATTGATTTCACGG 

m25 | L- 16314 CCACCCTTAACAGTACATAGTAC TDP62-57-52 140 
H-16409b (GTGCGGGATATTGATTTCACGG (5,5,25) 





*TDP — Touchdown protocol 


Table 2. Mitochondrial primers used for PCR amplification 








and aliquots of the extracted DNA. An alternative reaction mix, consisting of 100 mM dNTPs (25 
mM each dNTP) and 0.125—0.5 ul of Herculase( II Fusion DNA polymerase solution with the 1X 
Stratagene buffer, was used for the aDNA samples that were not successfully amplified using the 
AmpliTaq Gold reaction mix. The amplification cycling conditions were: 94°C for 5 minutes, then 
35 cycles of 94°C for 30 seconds, 52—63°C for 30 seconds, and 72°C for 30 seconds, followed by 
72°C for 10 minutes as a final extension, using a GeneAmp® PCR System 9700 Thermocycler 
(Applied Biosystems). 


Before proceeding with the cloning steps, excision and purification of the gel bands of specific 
sizes were executed using the Qiaquick® Gel Extraction Kit according to the manufacturer's 
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instructions. Cloning of the PCR products was carried out using the TOPO TA Cloning Kit and 
One Shot MAX Efficiency DH5a-T1 chemically competent E. coli cells (Invitrogen), following 
the standard protocols. Successful ligation and transformation were detected by the blue-white 
selection. Colony PCR was performed in a reaction mix of total volume 50p1, using the universal 
M13 primers, forward 5’-GTAAAACGACGGCCAG-3’ and reverse 5-CAGGAAACAGCTATGAC-3 
(Eurofins MWG Operon). The cycling protocol was 94°C for 10 minutes, and then 25 cycles of 
94°C for 1 minute, 55°C for 1 minute, 72°C for 1 minute, with a final extension at 72°C for 10 


minutes. 


For Sanger sequencing, ExoSAP-IT® (USB) was used for DNA purification. Cyclic sequencing 
for the colony PCR was carried out using a reaction mix containing 10 pmol of each M13 primer 
and 1X BigDye® Terminator v3.1 Cyclic Sequencing Kit (Applied Biosystems). The cycling 
conditions were 96°C for | minute, then 25 cycles of 96°C for 10 seconds, 50°C for 5 seconds and 
60°C for 4 minutes. Sodium dodecyl sulphate (Phytotechnology laboratories) was added to each 
PCR product, at a final concentration of 0.2% for further incubation at 98°C for 5 minutes, then 
25*C for 10 minutes and Centrisep columns were used for the filtration of the product prior to 
sequencing. A Genetic Analyzer (ABI Prism3130, Applied Biosystems) was used for the electrophoresis 
of purified amplicons, according to standard protocols. 


Diagnostic mutations in the mitochondrial sequences were determined using the majority rule. 
Sequence alterations that occurred with a frequency of more than 75% were defined as diagnostic 
mutations. Positions thought to be prone to post-mortem damage and artefacts were omitted from 
subsequent analyses (Gilbert et al. 2003, 32-47; Pusch and Bachmann 2004, 957—964; Gilbert et 
al. 2007, 1—10). 


Bioinformatic analysis: 


The Whit Athey's haplogroup predictor (Athey 2006, 34—39) [http://www.hprg.com/hapest5] was 
used to calculate percentages of the suggested Y-chromosomal haplotypes of the male mummies, 
(KV55/Akhenaten, Tutankhamun, Amenhotep III, Yuya, TT320-Cairo CG 61065) and to define 
the haplogroups of the ancient and modern control individuals. 


Matrilineal inheritance of mtDNA and characteristic single nucleotide polymorphisms (SNPs) of 
the royal family under investigation were determined using Lasergene (SeqMan and EditSeq) 
version 8.0 (DNAstar, Madison, Wisconsin). Mitochondrial DNA haplogroups were determined 
using the van Oven and Kayser phylogenetic tree (http://www.phylotree.org — mtDNA tree 2014; 
van Oven and Kayser 2009, E386—E394) and the Genebase mtDNA Haplogroup Reference Guide 
(www.genebase.com). 
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Results 


Y-chromosomal data: 


The Y-chromosomal data of the 7 male mummies are collectively shown in Table 3. Those subjects 
included 4 mummies from the Tutankhamun family group, and the father-son pair from another 
dynasty. One complete (KV55-Akhenaten) and 6 partial profiles (Amenhotep III, Tutankhamun, Yuya, 
TT320, Ramesses III, Unknown Man E) were obtained for the tested male subjects. On the other 
hand, all female mummies repeatedly tested negative. 


Using the Whit Athey haplogroup predictor, the most probable haplogroups obtained for the Tutankhamun 
family members were G2a for Yuya and R1b for KV55. Due to the incomplete profiles observed for 
Tutankhamun and Amenhotep III, analysis led to different probability figures, despite their concordant 
allele results. The haplogroup for these two mummies was thus predicted based on the full KV55 
data, particularly since the relationships were confirmed through a previous study (Hawass et al. 
2010, 638-647). The mummies that did not belong to the Tutankhamun lineage showed other haplotype 
probabilities (L and E1b1a). It should be noted that the identified R1b haplotype in Tutankhamun, 
KV55 and Amenhotep III was different from those detected in the kit control DNA 007 as well as in 
the staff members (1 and 3) as the profiles differed in several loci. 


Mitochondrial data: 


The mtDNA of 9 individuals was analyzed; 7 belonged to the Tutankhamun family (4 males, 3 
females) and 2 females who were not (Table 4). Using PCR primers specifically designed to 
amplify short mtDNA fragments, partial analysis of the HVRI (16056-16409) of the D loop 
control region was done. Diagnostic mutations in the sequences from the tested mummies were 
determined using the majority rule. The total number of sequences obtained to characterize the 
HVRI region was 631. Considering the alignment data from Thuya (generation I in the family), 
87 sequences were generated in total using 6 bone biopsies and 9 different primer combinations. 
Position 16224 was covered by 36 sequences from 10 different PCR products using 4 primer 
combinations. Of these 36 sequences, 32 showed a C and 4 showed a T residue. According to 
the majority rule and in the absence of contaminant sequences, this defines C to be the authentic 
base at mtDNA position 16224. Moreover, position 16311 was covered by 38 sequences in total 
and derived from 9 PCR products using 3 different combinations; 35 showed a C at this position, 
while 2 showed a T and 1 showed an A. 


Likewise for Tutankhamun (generation IV), 79 sequences were obtained in total using 4 bone 
biopsies and 8 different primer combinations. Position 16224 was covered by 51 sequences in total 
from 10 different PCR products using 5 different combinations. Out of 51 sequences, 49 showed 
a C, while 2 showed a T. Position 16311 was covered by 51 sequences in total from 10 different 
PCR products using 3 different primer combinations. Out of these 51 sequences, 50 showed a C 
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Mummy HVRI Diagnostic Mutation*,| Predicted Sequences | Min/ Max | Coverage 
tomb Frequency, Experiments!)  Haplogroups* per Alignment| Coverage 
Thuya 16224 (8996 - 4 exps) K 87 4/57 Full 
(KV46) 16311 (9296 - 3 exps) 

KV35EL (Tiye) 16224 (9096 - 5 exps) K 74 3/43 Full 

(KV35) 16311 (100% - 2 exps) 
KV35YL 16224 (8896 - 4 exps) K 72 6/55 Full 
(KV35) 16311 (10096 - 3 exps) 
KV55 (Akhenaten) 16224 (7896 - 4 exps) K 137 27/69 Full 
(KV55) 16311 (10096 - 3 exps) 
Tutankhamun 16224 (96% - 5 exps) K 79 2/65 Full 
KV62 16311 (9896 - 3 exps 
b- Control Mummies 
Mummy HVRI Diagnostic Mutation*,| Predicted Sequences | Min/Max | Coverage 
tomb Frequency, No. of Trials!) |Haplogroups*|per Alignment| Coverage 
Yuya 16224 (8096 - 2 exps) K 64 2/58 Full 
(KV46) 16311 (8996 - 3 exps) 

Amenhotep lll 16311 (10096 - 1 exp) H2b 80 17/47  |Partial 16118- 
(KV35) 16409 
Sitra-In None K excluded 18 - Partial 16118- 

KV60 16232 
Hatshepsut None K excluded 20 - Partial 16118- 
(KV60) 16232 
c- Staff Controls 
Staff Member HVRI Diagnostic Mutation Predicted 
Haplogroups* 
Staff 1 16189; 16192; 16223; 16292; 16355; 16519 L2a1 
Staff 2 16069; 16126; 16193; 16300; 16309 J1d 
Staff 3 16519 H (CRS) 
Staff 4 16304; 16519 H5 
Staff 5 16217; 16309 HV2 
Staff 6 16218; 16519 H20 
Staff 7 16069; 16126; 16145; 16187; 16217; 16261; 16300; 16311 J1b1 
Staff 8 16051; 16132; 16172; 16189; 16219; 16278 U6 
Staff 9 16069; 16126; 16214; 16231; 16305 J 
Staff 10 16093; 16224; 16304; 16311 Kila 
Staff 11 16129; 16309; 16318; 16362; 16397 U7 








*Alterations observed in the mitochondrial DNA with a frequency higher than 75% were designated as 
diagnostic mutations. 
*The phylogenetic tree (http://www. phylotree.org) and the Genebase mtDNA Haplogroup Reference Guide 
(www.genebase.com) were used to determine mtDNA haplotypes (Van Oven & Kayser 2009, E386-E394). 
t'Experiments" indicate the number of primer combinations used to generate differently sized PCR 
amplicons which harbour the diagnostic mutation. 
&Haplogroup determination for staff members was mainly based on the findings of Saunier et al. (2009, 
e97-e103) and Elmadawy et al. (2013, 338-341), in addition to general information from the van Oven & 


Kayser (2008, 


€386-e394) phylogenetic tree. 


[able 4. Characterization of mitochondrial DNA of 18^ Dynasty mummies 


while only 1 showed a T. Using the same rule, this also defines C to be the authentic base at positions 
16224 and 16311, respectively. 


The two diagnostic mutations observed in Tutankhamun's DNA sequences (Fig. 1) are indicative 
of the mitochondrial haplogroup K according to the van Oven and Kayser (2009, E386—E394) 
phylogenetic tree (http://www.phylotree.org) and the Genebase mtDNA Haplogroup Reference 
Guide (www.genebase.com). The haplogroup data for the other individuals is listed in Table 4. 
Based on the occurrence of the two diagnostic mutations 16224 and 16311, six members of the 
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Y:R1b 
Amenhotep III Tiye 
(KV35) (KV35EL) 
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mt:K 
Y: R1b 
Akhenaten ? 
(KV55) (KV35YL) 
mt: K 
Y: Rib 
Tutankhamun 


(KV62) 


Fig. 2. Family pedigree depicting mitochondrial and Y chromosome haplogroups 


Tutankhamun family group were characterized as haplogroup K, while one was characterized with 
haplogroup H2b. For the two KV60 reference mummies, a small PCR product encompassing 16118 
to 16232 was amplified and sequenced but it did not show the 16224 diagnostic mutation associated 
with haplogroup K. 


A pedigree representing the generated mitochondrial and Y-chromosome haplogroups of the investigated 
Tutankhamun family members is shown in Fig. 2. 


Discussion 

Implications from Y-chromosomal DNA analysis: 

Haplogroup R1b 

The proposed haplogroup for three of the investigated mummies, namely Tutankhamun, KV55 
(Akhenaten) and Amenhotep III, was R1b. The R1b haplogroup is defined by the mutation M343 and 


thought to originate in the western Asian region (Wood et al. 2005, 867—876; Cruciani et al. 2010, 
800—807; Myres et al. 2011, 95—101; Klyosov 2012, 87-105). R1b is not abundant in contemporary 
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Africa (Cruciani et al. 2010, 800—807), yet it is geographically widespread among the European 
populations making up to 50-60% of the genetic pool of the modern Europeans (Arredi et al. 2004, 
338-345; Cinnioglu et al. 2004, 127—148; Karafet et al. 2008, 830—838; Myres et al. 2011, 95-101; 
Klyosov 2012, 87—105). The introduction of this Y-chromosomal signature may be due to old 
introgression of lineages to the European gene pool from the Near East during the Neolithic spread 
of farming (Myres et al. 2011, 95-101). The phylogeographic analysis of the R1b haplogroup 
distribution provided strong support to the back-migration theory from Asia to Africa (Cruciani et 
al. 2010, 800—807). The R1b group was estimated to have an age of ca. 16,000 years before present 
(YBP) (Klyosov 2012, 87—105), and its suggested expansion time nearly approximates that of the 
G2a haplogroup (Cinnioglu et al. 2004, 127—148). Less than 7% of the contemporary Egyptians share 
the R1b Y-chromosomal fingerprint (Cruciani et al. 2010, 800—807; el-Sibai et al. 2009, 568—581). 


Haplogroup G2a 


The great-grandfather of Tutankhamun, Yuya, carries a Y-chromosomal signature that could be 
assigned to the haplogroup G2a. Haplogroup G is an F-affiliated clade (Luis et al. 2004, 532—544; 
Wood et al. 2005, 867—876), and it is defined by the mutation M201 (Cinnioglu et al. 2004, 127—148; 
Luis et al. 2004, 532—544; Wood et al. 2005, 867—876; Karafet et al. 2008, 830—838). This clade is 
not globally abundant, and its prevalence is mainly in the Middle East (highest in Druze), the 
Mediterranean basin and Caucasus Mountains where it exhibits its maximum frequency (Cinniošlu 
et al. 2004, 127-148; Karafet et al. 2008, 830—838; Balanovsky et al. 2011, 18255-18259; Lacan et 
al. 2011, 2905-2920). The pattern of this haplogroup distribution in the Caucasus suggests a Near 
Eastern origin (Cinnioglu et al. 2004, 127—148; Balanovsky et al. 2011, 18255-18259). The genetic 
share of the F-affiliated groups (G, H, I, J) is around 40% of the modern Egyptians, with G-M201 
representing approximately 9% of the population (Luis et al. 2004, 532—544). 


Haplogroup L 


One of the male control mummies, an unidentified 18^ Dynasty royal, previously thought to be 
Thutmose I (TT320-Cairo CG 61065), was proposed to belong to haplogroup L. This haplogroup is 
mainly defined by the M20 mutation, with additional markers to further confirm the clade definition 
(Sengupta et al. 2006, 202—221). It is found at low frequencies in some parts of the Middle East and 
Mediterranean Europe, and is mainly observed in southern, western and central Asia (highest in Indian 
subcontinent) (Kivisild et al. 2003, 313-332; Cordaux et al. 2004, 231—235; Sengupta et al. 2006, 
202-221; Thanseem et al. 2006, 42; Karafet et al. 2008, 830—838). In Y-chromosome studies on 
modern Egyptians, no haplogroup L findings were reported (Arredi et al. 2004, 338—345; Luis et al. 
2004, 532—544; el-Sibai et al. 2009, 568—581). 


Haplogroup E1b1a 


The father-son pair, Ramesses III and Unknown Man E/Pentawere, was found to have the haplogroup 
Elbla (Hawass et al. 2012, e8268) that shows its highest frequencies in modern populations from 
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West Africa (~80%) and Central Africa (—60%). It is less than 10% among North Africans and nearly 
absent in East Africa, which is where it was hypothesized to have originated (Trombetta et al. 2011, 
e16073). Alternatively, an earlier origin of around 20,000—30,000 YBP, was hypothesized in West 
Africa due to its very high diversity and frequency in that location (Rosa et al. 2007, 124). Although 
the frequency of haplogroup E among contemporary Egyptians is 39.595, that of the E1b1 subclade 
was only 1.4% (Cruciani et al. 2004, 1014—1022; Luis et al. 2004, 532—544; el-Sibai et al. 2009, 
568—581). 


Implications from mitochondrial DNA: 
Haplogroup K 


Haplogroup K most likely belonging to a Near Eastern lineage (Brandstátter et al. 2008, 191; Fernández 
et al. 2014, e1004401), is a subcluster of the macrohaplogroup U (Richards et al. 1998, 241—260; 
Kivisild et al. 1999, 1331—1334; Macaulay et al. 1999, 232—249; Maca-Meyer et al. 2001, 13). The 
age estimate of the K cluster is ~15,500—25,500 YBP (Richards et al. 2000, 1251-1276). The K 
haplogroup has been found to be present in 4.7% of contemporary Egyptian populations (Saunier et 
al. 2009, e97—e103). 


Haplogroup H 


The haplogroup of one of the studied mummies (Amenhotep III) was proposed as being H2b. 
Haplogroup H displays a sequence in the hypervariable region I (HVRI) identical to that of the 
Cambridge Reference Sequence (CRS) (Torroni et al. 1996, 1835-1850; Macaulay et al. 1999, 
232-249; Richards et al. 2000, 1251—1276; van Oven and Kayser 2009, e386—e394). It is believed 
that the HVRI 16311 transition is an emerging subclade of the H-CRS and named as H2b (Roostalu 
et al. 2007, 436—448; Brandstátter et al. 2008, 191). The age estimate of cluster H in Europe can be 
dated to 19,200—21,400 YBP; however, it is older in the Near East and can be dated back to 23,200— 
28,400 YBP. Haplogroup H is substantially prevalent in Europe as it is in the Near East but with a 
lower frequency (Torroni et al. 1996, 1835-1850; Richards et al. 2000, 1251—1276; Tambets et al. 
2000, 219—235; Malyarchuk et al. 2008, 1651—1658). This indicates that, although this cluster might 
have its roots in the Middle East, the significant evolution of this mitochondrial signature occurred 
in Europe (Richards et al. 1996, 185—203), probably in the late Upper Paleolithic expansions (14 500 
YBP) (Richards et al. 1998, 241—260; Richards et al. 2000, 1251—1276). The frequency of haplogroup 
H was, respectively, 3.6 and 5.94% among modern Egyptians (Saunier et al. 2009, e97—e103; Elmadawy 
et al. 2013, 338—341). 


Concluding remarks: 
This study reconfirms and complements the previously published family pedigree of the King 


Tutankhamun (Hawass et al. 2010, 638—647). It also presents the plausible mitochondrial and 
Y-chromosomal signatures of the family members. The data suggest geographical and temporal 
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matching ofthe origins of Y-chromosomal and the mitochondrial signatures in the family to possibly 
Near Eastern lineages, with time coverage believed to be between ~14,500 and 28,400 Y BP. A recent 
study by Schuenemann et al. (2017, 15694) has also shown the propensity of Near Eastern lineages 
among their ancient Egyptian sample coming from Faiyum, although these were from later ancient 
periods. A distinction between the geographical origins of the investigated late 18^ Dynasty family 
and the succeeding 20" Dynasty pharaohs (Ramesses III and Unknown Man E/Pentawere) was also 
shown. The latter was estimated to be of African origin, not a Middle or Near Eastern one. 


Due to its geographical location, Egypt could have over time created a melting pot of various gene 
pools. The genetic variation outside Africa reflects only a subset of the African gene pool (Underhill 
and Kivisild 2007, 539—564). Studies on today's Egyptian population samples have resulted in 
haplogroups that are predominantly European or west Eurasian (Stevanovitch et al. 2003, 23-39; 
Saunier et al. 2009, e97—e103; Elmadawy et al. 2013, 338—341; Fadhlaoui-Zid et al. 2013, e80293). 
It was also reported that the modern Egyptian haplogroups in a study were more closely associated 
with non-African populations than with the West Africans (Pagani et al. 2015, 986—991). This accords 
with our findings that the same trend was indicated through the analysis of 18^ Dynasty family members. 
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